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Abstract: We study the effect of shear on the cosmological backreaction in the context of 
matching voids and walls together using the exact inhomogeneous Lemaitre-Tolman-Bondi 
solution. Generalizing JCAP 1010 (2010) 021, we allow the size of the voids to be arbitrary 
and the densities of the voids and walls to vary in the range < Sl v < il w < 1. We derive 
the exact analytic result for the backreaction and consider its series expansion in powers 
of the ratio of the void size to the horizon size, ro/to- In addition, we deduce a very simple 
fitting formula for the backreaction with error less than 1% for voids up to sizes ro > in- We 
also construct an exact solution for a network of voids with different sizes and densities, 
leading to a non-zero relative variance of the expansion rate between the voids. While 
the leading order term of the backreaction for a single void- wall pair is of order (rn/io) 2 ; 
the relative variance between the different voids in the network is found to be of order 
(ro/to) 4 an d thus very small for voids of the observed size. Furthermore, we show that 
even for very large voids, the backreaction is suppressed by an order of magnitude relative 
to the estimate obtained by treating the walls and voids as disjoint Friedmann solutions. 
Whether the suppression of the backreaction due to the shear is just a consequence of the 
restrictions of the used exact models, or a generic feature, has to be addressed with more 
sophisticated solutions. 
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1. Introduction 

The non-commutativity of time evolution and spatial averaging in general relativity implies 
that the average expansion of a universe with structures does not evolve in time like the 
uniform Hubble expansion in a homogeneous Friedmann model [1—3]. This effect can be 
quantified by a backreaction term in generalized Friedmann equations and understood as 
a consequence of the nonlinearity of gravity [4] . Although the cosmological backreaction is 
conceptually well-understood, the complexity of the structure formation at the nonlinear 
level means its magnitude in the real universe is difficult to evaluate and is hence widely 
debated [5]: some studies have found that the backreaction is small [6—11], while other 
studies suggest it can have a significant effect on the cosmological dynamics and observa- 
tions [12-16], even to the extent of accounting for the observed cosmic acceleration entirely 
without additional effects [12-14]. As many of the current standard values for the cosmo- 
logical parameters - not just the cosmological constant - rely on the hypothesis of negligible 
backreaction, its evaluation has become a central issue in cosmology today [17, 18]. 

A widely used scheme to estimate the backreaction is to average the scalar parts of 
the Einstein equation on spacelike hyper surf aces, defined by constant proper time of the 
freely-falling dust particles. In this so-called Buchert approach, the backreaction term is 
given by the (positive) variance of the expansion rate minus the (positive) average shear [4]. 
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In studies that have estimated the backreaction to be significant [12-15], the shear on the 
boundaries between regions characterized by different expansion rates has been neglected 
or, equivalently, the boundary regions or matching conditions have been ignored. On the 
other hand, in perturbative studies that do not neglect the shear the backreaction has been 
found to be small [19-21]; however, see [22] for a discussion on the possible shortcomings 
of the perturbation theory in modelling the structure formation. Considering that on 
physical grounds, shear is expected to occur on the boundaries between regions of different 
expansion rates, evaluating the effect of the shear on the backreaction is evidently one of 
the key issues in the problem. 

This work is a continuation to our previous work [23], hereafter paper I, in which the 
effect of shear on the cosmological backreaction was studied in the context of matching 
voids (with £7 V = and Ht = 1) and walls (with S7 W = 1 and Ht = 2/3) together using the 
exact inhomogeneous Lemaitre-Tolman-Bondi or LTB solution. Whereas neglecting the 
exact matching can lead to significant backreaction, the main conclusion of paper I was 
that the shear arising from the exact matching suppresses the backreaction by the squared 
ratio of the void size to the horizon size, (ro/to) 2 , thus making it small for voids of the 
observed size ro/io ^ 1CP 2 [24]. Here we generalize the study of paper I by considering: 

1 . The backreaction as a function of the void and wall density parameters in the range 
< f2 v < f2 w < 1, thus relaxing the priors f2 v = and $7 W = 1 of paper I. 

2. The backreaction for voids of arbitrary size ro, thus relaxing the condition vq <C to 
assumed in paper I. 

3. A network of voids with different densities O v and radii ro, thus giving rise to relative 
variance of the expansion rate between the different void-wall pairs. 

Note that since the backreaction of the quasi-spherical Szekeres model reduces to the LTB 
model, the results of this work apply to the quasi-spherical Szekeres solution as well [25, 26]. 

Although we focus here on the effect of the structure formation on the average dy- 
namics of the universe, the ultimate goal of the research is to evaluate the total effect 
of the structure formation on the cosmological observations. In addition to the dynami- 
cal backreaction considered in this work, known effects of the structure formation include 
modifications on the propagation of light not directly determined by the volume-averaged 
expansion [8, 16, 27-41] and effects due to our non-average location in the universe [15, 42- 
44]. These effects may be important but we do not try to evaluate them in this work. 

The paper is organized as follows. The necessary background of the Buchert averaging 
method and the LTB solution are introduced in Sects. ^ and ||, respectively. Sect. |4] provides 
the definition of the void-wall LTB model, which we apply to study the effect of shear on 
the backreaction in Sect. ||: The general analytic expression of the backreaction for the 



void-wall LTB model is derived in Sect. 5.1 and its expansion as a power series and a simple 



but accurate fitting formula are considered in Sect. 5.2 . The result for the backreaction is 



then applied to a network of different voids in Sect. and to large voids in Sect. 5A. In 



Sect. 5.5, we discuss uncompensated voids. Finally, the results are summarized in Sect. 



and the conclusions are given in Sect. [?]. 
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2. Scalar averaging 



The spatial volume-average of a scalar S(x,t) is defined as 

L. S(x, t)e(x, t) 

where the volume element is determined by the determinant of the spatial part of the metric 
as e = y/det gijd^x and T> is the averaging domain or a region of the spatial sections. Unless 
otherwise noted, we take the averaging domain to be an origin-centered ball of coordinate 
radius R, i.e. V = M(R), and simply write (S) = (S)^jiy 

The volume expansion scalar 9 is defined in terms of the temporal change of the volume, 

9(x,t) = dtHe(x,t)) , (2.2) 

and is related to the generalized scale factor a(t) as 

<*>(*) = 3^, (2.3) 
a[t) 

where 

/ £ (x,t)N 1/3 



{2A) 



By taking the time derivative of the average (|2,l| ) and writing the expression in terms of 
the expansion scalar (2.2), we obtain the commutator of time evolution and averaging 

d t (S) = (d t S) + (S6)-(S)(8) , (2.5) 

which, when applied to the expansion scalar itself, yields 

d t (9) = (d t 9) + {(9-(9)) 2 \ . (2.6) 
>o 

The result ( |2.6| ) shows that the average expansion always decreases less (or increases more) 
in time than the average of the time derivative of the local expansion. 

By applying the averaging (^[l]) to the scalar parts of the Einstein equation in an 
irrotational dust universe, and using the results (2J3) and ([2.5|), we obtain the generalized 
Friedmann equations [4]: 

,S(t) 



a(t) 



-4nG(p)(t) + Q(t) (2.7) 



^|y) 2 = **GM) ~ ~ \m (2-8) 

dt(p)(t) = -3^<p)(i) , (2.9) 
a(t) 

where 1Z is the Ricci curvature scalar of the spatial sections and the backreaction Q is 
given by the variance of the expansion rate 9 minus the average of the shear scalar a 2 , 

Q(t) = I ((0 2 ) - (9) 2 ) - (a 2 ) . (2.10) 
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The backreaction ( 2.10| ) quantifies the difference of the time evolution of the averages 



relative to the uniform quantities in a homogeneous Friedmann dust universe. 

By partitioning the averaging domain D into a set of N mutually disjoint subregions 
T>i that satisfy (Ji=i = ^ the backreaction ( |2.1[1| ) can be written as 

N N N 

Qv = Y J fiQv I + -J2Y,f^(( e ^-(0)vf , (2.11) 

i=l i=l j=l 

where fi is the volume fraction of the subregion T>i, fi = Vol('Di)/Vol(P). The expression 



(2.11) makes it explicit that the backreaction is an intensive, rather than extensive, quantity 



and becomes useful in Sect. 5.3, where we consider a network of voids constructed by 



matching together different LTB solutions. 
3. LTB solution 

The exact spherically symmetric dust solution of general relativity was discovered by 
Lemaitre in 1933 [45] and is now commonly referred to as the LTB metric: 

ds 2 = _ dt 2 + l^M^dr 2 + A 2 (r, t)(d9 2 + sin 2 9^ 2 ) , (3.1) 
1 — k{r) 

where the prime stands for the radial derivative A'(r,t) = d r A(r,t), k(r) is related to the 
spatial Ricci curvature scalar as 

d r (A(r,t)k(r)) 
K ~ l A*(r,t)A>(r,t) ' [6 > 

and A(r, t) is determined by the Friedmann-like evolution equation, which, following the 
notation and parametrization introduced in Ref. [46], reads as 

1/2 



H(r,t)=H (r) 



°° (r) (iM)J + (l - n ° {r)) 



(3.3) 



where H(r,t) = dtA(r,t)/A(r,t) = A(r,t)/A(r,t), Hq(t) = H(r,to) and f^o( r ) are bound- 
ary condition functions specified on a spatial hypersurface t = to that determine the radial 
inhomogeneity profile, while the freedom to choose the function Ao(r) = A(r,to) corre- 
sponds to the scaling of the r-coordinate. The coordinate freedom is used here to set 

Mr) = r • (3.4) 

The curvature function k(r) in the metric Q3,l| ) is related to these by 

k(r)^H 2 (r)A 2 (r)(n (r)-l) (3.5) 

and the boundary condition function f2o( r ) is related to the physical matter density p(r,t) 
on the t = to hypersurface as 

n , 8ttG Jn (r) Po(r)d 3 x 



x 
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where po(r) = p(r,to) and d 3 x = r 2 sin6>drd6>d(/?. 

For general functions £lo(r) an d Hq(t) that are independent of each other, the LTB 
solution contains both decaying and growing inhomogeneities [47]. However, given the 
observed near isotropy of the CMB, the models with close to homogeneous early universe 
form perhaps the most relevant subcase of the LTB solutions. This corresponds to solutions 
where growing modes dominate so we only consider inhomogeneity profiles that obey the 
constraint [23]: 



Ho(r) = i 



1 - . (3.7) 



To calculate the backreaction fl2,lC| ) for the LTB solution in Sect. ||, we need the 
following quantities: the shear scalar 



ffM = ff v = 3 iM"^) 1 ' (3,8) 



the volume expansion scalar 



.,_ Mr 1 t) A'(r,t) _ d r (A\r,t)A(r,t)) 
1 ' 1 ~ A(r, t) + A'(r, t) ~ A 2 (r, t)A'(r, t) ' ( j 



and their expressions on the t = to hypersurface: 

a 2 (r,t ) = ^(rH' (r)) 2 , (3.10) 

6(r, t ) = 3H (r) + rH' (r) = ^d r {r 3 H (r)) . (3.11) 

For volume averaging, we also need the LTB volume element 

, s n , ,„ , A' (r,t)A 2 (r,t) sin 6 , 

e(x,t) = ^detgijdrdOdip = y ' y == drdOdp . (3.12) 

\/l-k(r) 

4. The void-wall LTB model 

We consider an LTB solution consisting of two different regions: a void with density param- 
eter fi v and a wall with density parameter f2 w , such that the boundary condition function 
( ]3,6| ) has the form 

tt (r) = + ( V 7 ^ " V^Mr - r )) * , (4.1) 

where G stands for the Heaviside step function and ro determines the size of the void. 
Because of the constraint (3/7), the density profile ( [P| ) implies the expansion profile 



+-i , 

H (r) = ^-(3- + (y^h - y^K)Q(r - r )j . (4.2) 
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To calculate the backreaction ( 2.10 ) for the void- wall model defined by Eqs. ( [4.1| ) and ( |4.2j ) 
in Sect. [|, we also need the first derivative of the expansion profile: 



K(r) = -foHV^ ~ V^KW - r ) , (4.3) 

where 5 stands for the Dirac delta function. 

Note that in the case O v = and O w = 1 the void-wall profile reduces to the step- 
function limit n — > oo of the LTB model considered in paper I. The reason to consider 
only the step-function profile is that, as already verified in paper I, the dependence of the 
backreaction on the sharpness of the transition between the void and the wall is weak. 
The step-function transition has also the advantage of allowing us to analytically perform 
calculations that would otherwise call for approximations or numerical methods. 



5. Backreaction in the void-wall LTB model 

In this section, we calculate the backreaction ( p.lO| ) for the void-wall LTB model defined 
by the density profile (|4.1| ) and the expansion profile (f4.2[) . Apart from the profile, the 
systematic study of the role of shear makes our approach different from the previous studies 
on the backreaction in the LTB model which have focused on finding profiles that exhibit 
acceleration of the average expansion [48-52] , on scale-dependence of the averages [53] or 
on general properties of the backreaction [54, 55]. 

When considering numerical values for the backreaction, we give the results in units 
of the backreaction obtained by taking the variance of two disjoint Friedmann models: the 
empty Milne solution for the void and the spatially flat Einstein-de-Sitter solution for the 
wall, yielding the result found in paper I: 



Qfrw = \tf > ( 5 - x ) 



where the regions are chosen to have equal volumes. These units are convenient as they 
directly address the role of shear by telling us how much the backreaction is suppressed 
relative to the Friedmann estimate which neglects the shear altogether. 

5.1 General expression 



With the help of Eqs. (3.8) and (3.9), the backreaction ( p. 10 ) of the LTB solution simplifies 
to 

oW _ 9 /^!M\ 4 a / Mr,t)A'(r,t) \ _ 2 / A(r,t) A'(r,t) \ 2 

^ t} ~ \A 2 (r,t) r \A(r,t)A'(r,t) 3 \ A(r,t) + A'(r,t) ' [ ' 



which, when evaluated on the t = to hypersurface using Eqs. ( |3.9| ) and ( |3.11 ), reduces to 



Q(fo) = 6 (<tf 2 } - (H ) 2 ) + 4 «ri^o> - (rflj) (H )) - \ (rH' ) 2 , (5.3) 



where Hq = Hq(t). 



- 6 - 



For the LTB solution, we write the volume average (2.1) as: 



" R r 2 dr 

, s . Jb^MH 7 ^) = Jo ^J^W) = [5] 
J m) e(r,t) f R r 2 dr ""wig 

7o V 1 -^^) 



(5.4) 



where we have defined the reduced dimensionless volume v, which, after the change of 
variables r = toy, reads as: 

rex y 2 dy 



V 



where e = ro/to, x = R/ro > 1, = 0(r — ro) and 



(5.5) 



r.((-)! = yQ v f( v () !V - v /'O v )(-) 



1 



r> v + (v^w - v^v)e . (5.6) 



To calculate the backreaction (|5.3|) , we thus need to evaluate v, [Hq], [Hq], [rH' ], and 
[rH' H ]. 

An integral appearing often in these quantities is 



which can be calculated analytically to yield 
1 



" y 2 dy 



(5.7) 



X(fi,v, a) 



2a 3 /2 



\fa.{v\J 1 + av 2 — \i\J 1 + afi 2 ) + In 



/ y/afi + a/1 + 



(5.8) 



For example the volume (5.5) can be written in terms of the function (|5.7|) as 



y 2 dy 



+ 



o yi + yMO) Je v/l + y 2 «(l 



y 2 dy 



X(0,e,a(0)) +X(e,ex,a(l)) , (5.9) 



where a is defined in Eq. (|5,6j). Using the definitions (5/7) and (|5.9|), [ffo] and [-ffg] can be 
integrated to yield 



[fib] =tU[l 



2 + 1 



[fio] = *o 4 1 



2b + 1 



(5.10) 
(5.11) 



which combine to give 



<fio 2 > " (Ho) 



\2 



1 (VTl^ - y/SQ 



-2q2i 



(5.12) 
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Similarly, using Eq. (|4.3| ), we have for the derivative terms 



[rH' ] 
[rH Q H' } 



-^(v/JU-VJT 



„,3 



1 1 5 (t y - r )dy 
o s/T+yMQ) 



(5.13) 



/a 



£Z „,3 



V ' 5 (toy - r )dy 
o s/T+yM@) 



+ 



ex y 3 @(t y - r )5(t y - r )dy 

Vi + y 2 a(e) 



(5.14) 



The integrals in Eqs. ( |5.13 ) and ( j5.14| ) can be recognized as Stieltjes integrals: using 

d@(t y - r 



dy 



toS(t y - r ) 



they can be written as the following ordinary integrals 

dtf 



A 
B 



t- 

Jo Jl 



^Jl + e 2 a(d) 
1 $d$ 
o y/l + e 2 a(#) 



so that Eqs. ( |5.13D and (|5.14|) become 



\rH'} 



[rH H' ] 



3 



to 



-e 



Vi Y )A 

3 ( ■ - \/Q 



/a 



A — -l v'O 



(5.15) 

(5.16) 
(5.17) 

(5.18) 
ft v )/n. (5.19) 



By recalling the definition (5.4) and substituting the expressions ( 5.1CQ , ( |5.12 ), ( 5.18j ) 
and ( |5.19|) in Eq. ( |5 . 3|) , we finally obtain: 



Q = t 



-X Xi + - [ (Jo + Xi) B - Xa] e 3 - — A 2 e 6 



(5.20) 



where the integrals ( f>. 16 ) and ( |5.17 ) can be calculated numerically or as an expansion in 
powers of e or e _1 . 

5.2 Power expansions and a fitting formula 

As the cosmological observations suggest that voids in the cosmic matter distribution satisfy 
ro <C to [24, 56-59], one of the most useful methods to investigate the backreaction is the 
series expansion in powers of e = ro/io- By expanding the exact result ( 5.20| ) in powers of 
s to second order, we obtain: 

Q = -^tf (g) 3 (nj/ 3 - ^ /2 ) 2 (54 - nl - ml - i5n w n v - 6^y 2 + 



2on^ 2 nl! 2 + 45fi w ^y 2 + 9on v nH 2 - sonl^nl! 2 + 36ft 3/2 + 



+ 9ft 3 / 2 - 20ft w - 48ft v - 30fti /2 - 36^y 2 ) f 7 ^) + O 



(5.21) 
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while for the values V = and f2 w = 1 the fourth order expansion reads as 

2 



i? 



135 \*o 



2677 



3359 /r \ 3 



102060 437400 \RJ 



where the first term in the expansion agrees with the result given by Eq. (5.11) in paper I. 
For the values f2 v = 0, $7 W = 1 and (ro/R) 3 = 1/2, the sixth order expansion of Eq. ( |5.2C| ) 
yields: 



135 V t 



137107 



12247200 V t 



ro 



33336241 



3940536600 \t 



r_oY_ (r_o 



to 



(5.23) 



By inspecting the power series ( 5.23| ), we see that a very simple but accurate fitting 
formula for the backreaction in terms of elementary functions is given by 



Q 



l-2 



135 



1 



4 V^o 



(5.24) 



Testing the fitting formula ( 5.24j ) numerically shows that it is very accurate up to horizon 
sized voids ro ~ to (for r o < *o we have | error | < 1%) and an excellent approximation even 
beyond. To illustrate this, we have plotted the fitting formula fl5.24| ) against the exact 
backreaction ( |5,20D and the leading order term of the expansion (|5.23| ) in Figure |l[ The 
figure also shows that the mere leading order term gives an accurate approximation for the 
backreaction even up to voids of size ro ~ to/3. 




Figure 1: The backreaction as a function of ro/to for the profile (^1) with tt v = 0, ft w = 1 and 
(ro/R) 3 = 1/2 calculated using the exact result ( |5.20| ) (red solid curve), the leading order term of 
the expansion ([[2|) (blue dash dotted curve) and the fitting formula d5.24| ) (black dashed curve). 
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By solving for the density parameters fi v and f2 w that yield the maximum value for 
the backreaction of small voids ( |5.21 ), we obtain Q v = and f2 w = 0.7. This can be seen 
in Figure |2[ where the backreaction is plotted as a function of O w for three values of ro / to 
and fixed O v = and (ro/i?) 3 = 1/2. The figure also shows that the backreaction is only 
slightly larger for Q w = 0.7 than for $7 W = 1. 




Figure 2: The backreaction as a function of f2 w for the profile ( j4.l|) with ro/io = 0.01 (blue dash 
dotted curve), ro/io = 0.02 (red dashed curve) and ro/io = 0.03 (black solid curve). All profiles in 
this figure have O v = and (ro/R) 3 = 1/2. 



5.3 Network of voids 

Using LTB solutions with the profile defined by Eqs. (4.1) and (4.2), we can join together 



many void- wall pairs to construct a model for a network of voids with different f2 v and ro, 
presuming the following matching conditions are met: the different void- wall pairs must 
have the same wall-density f2 w and the same age of the universe tQ. This provides an exact 
solution, because, outside the void at R > ro, the void- wall LTB profiles are identical to 
the homogeneous Friedmann dust solution with f2 w as the density parameter and to as the 
age of the universe and thus naturally match together. 



For a network of voids, the total backreaction Q tot is given by the formula (2.11): the 
volume-weighted average of the backreactions Qi of the individual void-wall pairs i plus 
a sum over the relative variances of the expansion rates between all the different pairs. 
The relative variance term makes the configuration particularly interesting as it seems to 
offer a way to increase the variance of the expansion rate without having to introduce any 
counterbalancing extra shear. 

Let us determine an upper limit to the relative variance term for a network of voids. 
As the relative variance is given by the double sum term in the expression ( |2.11| ), we need 
to calculate the volume-average expansion (8) for a void-wall LTB model as a function of 
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the parameters f2 v , f2 w , ro, to an d -R: 

(0) = (3 - y^)t 1 + ^Uo 1 (^) 3 (^-y^)(l08 + 12^_ 6 ^3/2y^ + 



+ 18 - 90\/n^ + 45 y^^w + 45 V / ^^v - 80fi w - 80\/^v^w + 640 v + 



63a 3 / 2 - 6\/aa 3 / 2 - 6n w n v + 450?/ 2 - 6fl 2 



(5.25) 



Given the matching condition 0$ = 0,^' between different void- wall LTB regions i 
and j, the leading order terms of the volume expansion ( [5,25 ) cancel in the difference 
(9}i — (0)j. Consequently, the relative variance term in Eq. (|2.11 ) gets the maximum value 



>(j) 



when £1 



and ft 



U) 



1, yielding the upper limit: 



«*>« - (ow 2 < 



23_.-t fro 
540 U 



Therefore we have 

N N 



|EEm(^-^) 2 ^EE 



N N 1 son / \ 4 

1 529 2 ( r o 

En 



i=l j = l 



or 



3 N 2 2916OO V tn 

t=l j=l 

5EE/i/i((^-<^) a <io-v(?) 4 . 

1=1 1 = 1 



V TV 



(5.26) 



(5.27) 



(5.28) 



where ro is the radius of the largest void in the network. Eq. ( 5.2Sj ) shows that the relative 
variance term is at most of the order (ro/io) 4 so the backreaction for a network of LTB voids 
with profile defined by Eqs. Q4.1| ) and (|4.2| ) is essentially no greater than the backreaction 
for a single void- wall pair. 

5.4 Large voids 

Let us then study the backreaction of large voids, i.e. voids that do not satisfy ro <C to. 
For voids of the horizon size ro ~ to, the exact backreaction (5.20) has to be evaluated 
numerically or using the fitting formula ( 5.24| ) but for superhorizon voids with ro 3> to, we 
can consider the limit ro/io ~~ ► 00 of Eq. (5.20), which, after some manipulations, takes the 
form: 

2 (x 2 - 1) 8 
27 



Qc 



t, 



" U/^ + v^ 2 " 1 )) 2 
y/a^ + ^/o^{x 2 — 1 



+ 



9 



3 ^a v a v 
1 tfdtf 8 (x 2 



d-d 



1 



1 dtf 



+ 



(5.29) 



/o yja{d) 9 ^/a^ 7 

where a($) is defined in Eq. (|5.6j ), a v = a(0), (1) and the integrals can be evaluated 

in terms of elementary functions: 

'3V2(tf-i) N ^ 



dt? 



3^/2 



o y^MWj 4(vn; 



/a 



■ arctan 



4 y 7 !^ 2 " 



(5.30) 
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+ 



3\/2, 



3 ( arcsin v/f2 v — arcsin \l Q 



(3 — y fi v ) arctan 



^ 4 vr^ 7 ^ 



(5.31) 



Eq. ( 5.29| ) can be used to solve the values of the density parameters f2 v and $7 W that 
maximize the backreaction for the superhorizon voids: numerically, we find that the max- 
imum is located at S7 V = and fl w = 0.625 and has the value 



Qr 



0.038tn = 0.23Q F rw 



(5.32) 



which can also be read off from Figure |3[ The location of the maximum differs only slightly 
from the point Q v = , il w = 0.7 that gives the maximum backreaction for small ro < 0.3 to 
voids, although the maximum value of the backreaction relative to the value at fi w = 1 is in 
this case noticeably larger (cf. Sect. |5.2| ). However, even with the superhorizon sized voids, 
the maximum backreaction ( |5.32 ) is still only ~ 20% of the FRW value (5.1), implying 



that in the void- wall LTB models with the profile defined by Eqs. ([O]) and ( |4.2| ) the shear 
suppresses the backreaction in all cases at least by about an order of magnitude. 

In Figure |], we have plotted the backreaction as a function of the size of the void ro /to 
for different values of the wall density £l w but keeping the void density and the volume 
fraction fixed: Sl v = and (r /R) 3 = 1/2. The fi gure illustrates how the rapid growth of 
the backreaction as a function of the void size ro stops once the horizon size to is exceeded. 



0.20 



0.15 




Figure 3: The backreaction as a function of J7 W for the profile ( J4. l| ) with ro/io = 1 (blue dash 
dotted curve), r /t Q = 3 (red dashed curve) and r /to — > 00 (black solid curve). All profiles in this 
figure have fi v = and (r /R) 3 = 1/2. 
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Figure 4: The backreaction as a function of ro/to for the profile (4.1) with f2 w = 1 (black solid 
curve), O w = 0.3 (blue dash dotted curve) and fi w = 0.62 (red dashed curve). All profiles in this 
figure have fi v = and (r /i?) 3 = 1/2. 



5.5 Uncompensated voids 



As the density function Oo(r) is determined by the integral of the physical matter density 
( |3.6| ), the profile (4.1) implies an overdense shell in po(r) between the void and the wall (see 
Sect. 5.4 of paper I for a more detailed explanation). To avoid the overdense or collapsing 
shell, we consider the profile defined by 



9(r,to)=to l (3-Q(r-r )) 



(5.33) 



which, by virtue of Eqs. ( |3.7| ) and ( p.ll ), implies 



"o(0 = (i-(?) 3 ) 2 e 



(r - r ) 



(5.34) 



and 



H (r) 



ro\3 
r / 



G(r - r ) 



(5.35) 



In paper I, we found that the difference in the backreaction between the compensated 
and the uncompensated void was 0(10%) for small voids. In Figure ||, we have plotted the 
backreaction for the profile defined by Eqs. ( 5.34 ) and ( 5.35| ) to show that the difference 
remains small even for large voids. This proves that the suppression of the backreaction is 
not due to the compensating overdensity but truly an effect of the shear. 
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0.2 0.4 0.6 0.8 



Figure 5: The backreaction as a function of ro/to for the compensated void- wall profile ( [4.1[ ) with 
Q v = and fl w = 1 (red dashed curve) and for the uncompensated void-wall profile ( |5.33| ) (black 
solid curve). Both profiles in this figure have (ro/R) 3 = 1/2. 



6. Summary of the results 

We have generalized the analytic approach of Ref. [23] or paper I in studying the role of 
shear in the cosmological backreaction problem. As in paper I, we constructed a void-wall 
model from the exact inhomogeneous LTB dust solution, but instead of assuming the voids 
to be small ro <C to and fixing the densities in the voids and walls to O v = and O w = 1 , we 
let the void size ro be arbitrary and the densities vary within the range < O v < f2 w < 1. 
Moreover, we constructed an exact solution for a network of voids with different densities 
Sl v and radii ro to study how the emerging relative variance of the expansion rate between 
the different void-wall pairs affects the backreaction. 

We calculated the exact analytic result for the backreaction of a void-wall pair with 
profile defined by Eqs. (4.1) and (|4.2j ) in Sect. |5.1[ To study its general behavior for 



subhorizon voids, we considered the series expansion in powers of ro/io in Sect. 5.2. The 
leading order term of order (ro/to) 2 in the expansion was verified to agree with the result 
obtained in paper I for the case O v = and £l w = 1. From the leading order term we found 
that the values Q v = and £l w = 0.7 yield the maximum value for the backreaction for 
subhorizon voids of size ro < to/3. However, the increase in the backreaction was showed to 
be only ~ 10% relative to the case where Sl v = and fi w = 1. Furthermore, by inspecting 
the series expansion with the values fi v = 0, fi w = 1 and (ro/R) 3 = 1/2, we were able 
to deduce the simple fitting formula ( [5.24 ) for the backreaction in terms of elementary 



functions, with error less than 1% up to horizon sized voids. 



In Sect. 5.3, we considered a network of subhorizon voids and showed that the relative 



variance of the expansion rate between the different voids is at most of order (ro/io) 4 - The 
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total backreaction for the network is then essentially given just by the volume- weighted 
average of the backreactions of the individual void-wall pairs and thus remains of order 
(ro/io) 2 ° r small for voids of the observed size. 



We investigated the behavior of the backreaction for large voids in Sect. 5.4, in par 



ticular super horizon ro S> to voids by considering the limit ro/to ~~ * 00 • We found that the 
backreaction approaches a constant value as the size of a super horizon void is increased, 
the limiting value depending on the parameters f2 v , O w and R/r^. The ultimate maximum 
value for the backreaction was found to be Q m ax = 

0.038i Q 2 at O v = and Q w = 0.62 in 
the limit ro/to ~~ ► °°> keeping the void-to-wall volume ratio fixed to (ro/i?) 3 = 1/2. This 
is still only ~ 20% of the maximum backreaction obtained when the shear is neglected by 
treating the voids and walls as separate Friedmann solutions. Furthermore, as opposed to 
the ~ 10% increase in the case of small voids, the backreaction for large voids is a few 
times (depending on the ratio r^/to) greater at the maximum point Q v = and J7 W = 0.62 
than at fi v = and fi w = 1. Finally, in Sect. [5l] we demonstrated that the suppression of 
the backreaction is not due to an overdense or collapsing region between the void and the 
wall but truly an effect of the shear. 



7. Conclusion 

By considering exact solutions of the Einstein equation consisting of one or more LTB 
solutions, we have pinpointed the issue of small versus large cosmological backreaction to 
the question of matching conditions: while the variance of the expansion rate alone can 
induce significant backreaction, the shear arising from matching together the regions with 
different expansion rates seems to bring down the backreaction by at least five orders of 
magnitude for voids of the observed size. The crucial question is whether the suppression 
of the backreaction due to the shear is a general property of all realistic cosmological 
solutions of general relativity or just a special property of the matching in the considered 
particular solutions. This issue has to be addressed with solutions more sophisticated than 
the LTB-based models employed here. 
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